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Summary 


A  conformal  transformation  formula  using  Riemann-Stieltjes 
integrals  is  derived  for  use  with  problems  involving  the 
interaction  between  a  given  finite-sized  geometry  and  a  known 
far  field.  The  derivative  of  this  trans formation  is  non-singular 
in  the  domain  considered  and  tends  to  one  at  infinity.  A  formula 
is  derived  for  transformation  from  the  unit  circle  to  the  exterior 
of  an  arbitrarily  given  continuous  curve  with  bounded  variation. 

A  special  case  of  the  transformation  is  very  similar  to  that  of 
Schwarz-Christoffel.  Application  to  the  generation  of  aerofoils 
gives  some  fairly  flexible  formulae  of  the  finite  product  type. ^ 
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INTRODUCTION 


A  number  of  problems  involving  the  interaction  between  a 
given  finite-sized  geometry  and  a  far  field  such  as  potential  flows 
about  aerofoil  sections  (1],  [2],  {3],  scattering  of  two  dimensional 
long  waves  (4],  {5]  and  stress  concentration  around  two  dimensional 
holes  [6] ,  [7]  can  all  be  neatly  solved  using  external  conformal 
mappings.  Their  solutions  are  derived  from  those  for  the  simple 
basic  circular  geometry.  Therefore  the  common  problem  is  to  find  the 
conformal  transformation  from  the  unit  circle  | 5 |  ■  1  to  a  given 
closed  curve  in  the  z  plane  [8] .  This  is  equivalent  to  determining 
all  the  equipotentials  and  force  lines  external  to  the  given  closed 
curve.  Solution  of  the  potential  problem  using  the  greatest  lower 
bound  of  superharmonic  functions  19] ,  [10]  is  not  a  practically 
attractive  method  unless  the  field  is  discretized  and  handled  in  a 
Finite  Difference  scheme.  This  leaves  conformal  mapping  as  the  only 
direct  method.  However  the  required  conformal  transformation  must 
satisfy  the  following  conditions. 

(a)  The  mapping  must  be  one-to-one  onto  between  the 

regions  under  consideration.  _  (li 

(b)  The  mapping  must  be  continuous  and  conformal 

between  the  two  regions.  (1] 

(c)  The  difference  (z  -  5)  must  vanish  to  order 

—v 

C  at  infinity,  where  v  is  a  positive  number.  (It 

These  conditions  are  illustrated  by  examining  the  following 
familiar  transformations: 


The  Joukowski  transformation 

z  *  (z  -  m) (?  -  n) ((  -  (m  +  nJ/2)”1  (2] 

where  m  ■  1,  |n|  <1  satisfies  all  three  requirements  and  is  widely 
used  in  theoretical  aerodynamics.  Likewise  a  slight  change  in  the 
roots  m,  n  turns  the  curve  in  the  z  plane  into  an  ellipse  which  is 
also  widely  used  in  elastostatic  and  elastodynamics. 

Finite  trailing  edge  aerofoils  and  lenticular  shapes 
bounded  by  two  circular  arcs  are  generated  by  the  K£rm£n  -  Trefftz 
formula  ([11],  page  78) 

z  ♦  v  .  , 

z-v  '5-1'  ^ 

where  v  is  a  positive  number  0<vS2.  Again  the  three  requirements  are 


satisfied.  Two  cuts  are  required,  between  t  *  -1 ,  £  -  +1  and  between 
z  ■  -v,  z  ■  +v,  for  the  £  plane  and  z  plane  respectively.  The  sucess 
of  these  two  transformations  is  in  part  owing  to  the  containment  of 
all  singularities  of  dz/d(  inside  the  circles  in  the  C  plane. 

The  Schwarz  -  Chris toff el  transformation  S(o),  ([10],  page 
236)  defined  by 

dS  ^  _  2  ft  ^ 

•T-  -M Il(o-o.  )  Dk,  £  6  -  1 ,  M  a  complex  constant,  (4) 

do  .  k  .  ,  K 
k*l  k“l 

is  general  enough  to  generate  a  number  of  polygonal  curves  but  it 
does  not  satisfy  the  condition  (lc) ,  and  therefore  is  not  appropriate  to 
the  present  problems.  It  should  be  noted  that  if  the  o^'s  are  placed 
at  the  vertices  of  a  regular  n-aided  polygon,  all  0^  are  equal  to 
1/n,  then  we  get  a  transformation  from  the  circle  passing  through  n 
points  o^'s  in  the  o  plane  to  a  regular  n-sided  polygon  in  the  z 
plane.  As  the  number  n  tends  to  infinity  the  Schwarz  -  Christoff el 
transformation  becomes 

dS  m 

do  “  o*  ' 

which  is  a  mapping  from  the  outside  of  one  circle  to  the  inside 
of  another.  A  modification  of  this  transformation  will  be  considered 
in  Section  5. 


Another  transformation  method  uses  the  infinite  series 
*<C>  -  c  +  a^"1  +  a2?'2  +  ...  (5) 

and  can  generate  any  shape  as  the  series  is  a  Laurent  expansion  for 
the  complex  function  z(C)  which  is  analytic  outside  the  unit  circle 
|c|  ■  1.  However  it  requires  quite  a  number  of  terms  even  for  the 
simplest  aerofoil  sections  and  yields  no  information  whatsoever  about 
the  singularities  or  one-to-one  onto  property.  It  is  noted  that  Glauert 
[11]  has  applied  this  method  to  the  problems  of  aerofoils.  A  well 
developed  theory  on  the  determination  of  these  coefficients  for 
aerofoils  is  due  to  Theodorsen  and  can  be  found  in  [12] . 

Here  we  propose  a  transformation  of  the  kind  given  below 
by  equation  (6)  or  (10)  subjected  to  the  condition  (7) .  This 
transformation  will  be  shown  to  be  indeed  the  general  form  for  any 
conformal  mapping  satisfying  the  three  conditions  (1) .  A  formula 


[33 


is  then  derived  for  transformation  from  the  unit  circle  to  an  arbitrary 
given  shape.  A  test  case  is  given  for  transformations  of  the  finite 
product  type.  One  special  case  of  this  transformation  is  very  similar 
to  that  of  Schwarz  and  Christoff el.  Application  to  the  problem  of 
aerofoil  sections  gives  some  fairly  flexible  formulae  of  the  finite 
product  type. 


2.  ANALYSIS  OF  THE  NEW  TRANSFORMATION 

Let  t  be  a  real  variable  (OStSl) ,  a(t)  be  a  periodic  complex 
valued  function  of  t  with  bounded  variation  and  $(t)  (|$(t)|i  1) 
be  a  periodic  continuous  con^plex -valued  function  of  t.  He  have  the 
function 

*(C>  -  /^(t)do(t)  +  /*♦< t)dt  +  exp[/^log (c— ♦(t }  )da (t)  +  J*log(c-*(t)  »t] 

(6) 

defined  for  all  |(|  >  1  as  all  the  conditions  for  the  Riemann  -  Stieltjes 
integrals  are  satisfied.  The  function  log  (£-+(t))  is  chosen  such 
that  it  is  continuous  for  OStSl  with  each  given  value  of  £.  Furthermore, 
if  the  images  of  the  circles  |?|  ■  1  +  e,  e  >  0,  e  very  small,  are  non¬ 
self-intersecting  closed  paths  and  also  if  it* is  assumed  that 

J*  d  (a(t)  +  t)  *  0  for  all  finite  c,  |c|  >  1  (7) 

then  all  the  requirements  (1)  are  satisfied.  This  significance  of  the 
assumption  (7)  is  discussed  in  detail  below,  in  Section  4. 

First  z(£)  is  continuous,  its  derivative  is 

gjr  -  exp  [/*log(c-$(t))d(a(t)  +  t)]/*  d(a(t)  +  t) ,  (8) 

and  z(?)  is  thus  an  analytic  function.  The  exponential  factor  is 
always  non-zero  and  the  last  integral  has  been  assumed  in  equation 
(7)  to  be  non-zero.  Hence  condition  (lb)  is  satisfied. 


z  - 


For  large  value  of  5  we  have  the  expansion 
/^(t)d(a(t)  +  t)  ■  c  exp|/*log(l  -  -$-£!■  )d(a(t)  +  t)J 

-  ;  exp  {  -  !  1  /V(t)d(a(t)  +  t)  \  , 

^  k*l  ,  k  °  J 


(9) 


using  Arzela's  theorem.  After  an  expansion  for  the  exponential 
function  we  get  the  condition  (lc)  satisfied,  it  is  then  easy 
to  show  that  for  large  enough  value  of  |;|  the  mapping  z(t)  is  one- 


to-one 


The  remaining  part  of  the  proof  is  to  prove  that  the 
mapping  z(C)  is  one-to-one  onto  between  the  exterior  of  the  unit 
circle  in  the  ;  plane  and  the  exterior  of  its  image  in  the  z  plane. 

Now  let  c  describes  a  circle  of  very  large  radius  centered 
on  the  origin  of  the  (  pl^ne.  The  image  z(?)  describes  a  closed 
path  in  the  z  plane  which  is  very  close  to  this  large  circle  if  we 
superimpose  the  two  planes.  Let  the  large  circle  in  the  c  plane 
expand  continuously  as  |(|  tends  to  infinity.  The  closed  path  in  the 
z  plane  then  also  expands  continuously  and  sweeps  through  every  point 
exterior  to  its  initial  state.  Hence  outside  some  circle  |(|  *  M, 
the  mapping  z(t)  is  one-to-one  onto. 

Since  the  mapping  is  conformal  with  dz/d(  non-zero  everywhere 
we  can  unambiguouslydefine  the  "interior"  and  "exterior"  of  each 
closed  path  which  is  the  image  of  some  circle  centred  on  the  origin 
in  the  (  plane.  Now  draw  a  nest  of  concentric  circles  in  the  (  plane 
between  |c|  ■  1  +  e  and  |(|  «  M.  The  image  of  this  nest  is  another  nest 

of  closed  paths  in  the  £  plane.  Each  closed  path  of  this  nest  is  wholly  on  the 

interior  or  exterior  of  every  other  of  the  nest,  and  it  can  have  closed 
paths  as  close  to  it  as  we  like  on  both  its  exterior  or  interior,  we  now  draw 

a  circle  |(|  *  M  in  the  z  plane  and  let  it  shrink  continuously  to  ■  1  +  e. 

Its  image  closed  path  also  shrinks  continuously  from  its  initial  state 
to  its  final  state.  No  region  on  the  z  plane  bounded  by  these  two 
closed  curves  can  be  left  outside  this  sweep  by  consideration  of  the 
properties  of  the  nest  of  image  closed  paths  in  the  z  plane. 

The  mapping  z(()  is  also  one-to-one  in  the  range  M2|c|kl  +  e. 

For  points  of  different  radii  in  the  C  plane  we  have  different  closed 
paths  in  the  z  plane.  The  image  of  each  circle  centered  on  the  origin 
of  the  (  plane  can  not  fold  back  on  itself  either,  as  dz/d;  is 
everywhere  non-zero.  Therefore  no  two  points  in  the  (  plane  between 
the  circles  | £ |  ■  1  +  e  and  | ( |  ■  M  can  have  the  same  image . 

This  completes  the  proof. 

The  function  z({)  for  |(|  ■  1  is  defined  to  be  the  limit 
as  jc|-*l,  |e|>l  of  z(c).  Note  that  due  to  the  one-to-one  property 
of  the  mapping,  |z(0|  is  bounded  as  |c|  tends  to  unity. 

When  $(t)  assumes  the  value  exp(i2irt),  we  have  a  special 
case  of  equation  (6) ,  which  is 

z(0  -  J*  el2wt  do(t)  +  exp  [J1  log(c  -  ei2irt)  fda(t)  +  dtj] .  (10) 

9  O  'O 


15] 


With  an  appropriate  cut  in  the  C  plane,  we  can  show  that  equation 
(10)  is  indeed  the  general  fora  for  any  conformal  transformation 
satisfying  conditions  (1),  as  in  the  following; 

Let  W(c)  be  a  given  conformal  transformation  satisfying 
conditions  (1)  and  W(el2ir^)  is  a  function  of  t  with  bounded  variation 


Choose  a  point  WQ  strictly  interior  to  the  image  of  the 
unit  circle  |(|  >1.  Consider  the  contour  integral 

f  log  (W(U  -  Wc)^2_  -  o 

with  the  cut  on  the  C  plane  following  the  unit  circle  u  then  radiating 

from  the  point  +1  to  +•.  The  curve  y  is  along  the  unit  circle,  the 

radial  cut,  the  infinite  circle,  and  back  to  the  unit  circle  along  the 

radial  cut.  In  this  way  the  function  log  (W(C)  -  W  )  is  defined  and 

o 

single-valued  for  all  points  C  not  on  the  cut.  The  components  of  this 
contour  integral  give 


-2iri  log  (W(C)  -  Wc)  -  ^log  tW(o)  -  Wc]  — +  Jj  2iri 


-  \  -  V  rh-  • 

'*ere  a  c^rcle  of  radius  R  centered  on  the  origin,  R  tending 

to  infinity. 


The  first  integral  on  the  right  hand  side  is  then  integrated 
by  parts.  The  third  integral  has  its  integrand  expanded  using  the 
condition  (lc) .  The  result  is 


log  (W(c)  -  W  )  -  J1  log  <c  -  ei2irt)d 
C  t-0 


log  (W(ei2irt)  -  w  ) 
c 

2iri 


The  square  bracket  is  a  continuous  complex- valued  function 

^  2  irt 

of  t.  Due  to  our  assumption  that  W(e  )  is  of  bounded  variation 
the  square  bracket  in  the  last  integral  also  has  bounded  variation. 
Therefore 


W(?)  -  W  «  exp (J1  log  (c  -  ei2wt 
t-0 


)d(6(t)  ♦  t)) 


where  B(t)  is  a  continuous  periodic,  complex  valued  function  of  t, 
having  bounded  variation.  (The  function  B(t)  may  not  have  bounded 
variation  if  V?c  is  on  the  image  of  the  unit  circle) . 

Routine  integration  shows  that 

Co  -12,t  a 

satisfying  the  requirement 

W c  *  /  e  dB(t) 

0 

of  equation  (6) . 

Hence  formula  (10)  is  indeed  the  general  form  for  any 
conformal  transformation  z(()  satisfying  the  conditions  (1)  subjected 
to  the  conditions  (7)  and  that  W(el27Tt)  is  of  bounded  variation. 

Section  3  will  show  that  equation  (6)  is  only  a  special 
case  of  equation  (10) ,  both  subjected  to  the  condition  (7) .  Henceforth 
we  take  (10)  as  the  standard  form  for  our  transformation. 

a 

The  form  of  equation  (10)  can  be  simplified  even  further 
when  the  imaginary  part  of  the  function  a(t)  is  continuous.  Let 
us  denote  Xm(a(t))by  b  (e1  j  and  consider  the  integral 

r1  ,  i2irt.  „  .  i2irt.  f1  .  ,  i2itt,_  ,  .  i2irt. 

/t-0  109(5  "  *  )db(e  )  "  "  Jf o  b(e  )d  109(5  “  *  > 

which  is  the  imaginary  part  of  the  second  integral  in  equation  (10) 

Since  b(e  )  is  continuous,  the  internal  boundary  value  problem  for 
the  unit  circle  has  a  solution  b(£  +  in)  which  is  harmonic  in  the  (€,  n) 
plane  and  takes  the  given  boundary  values  on  the  unit  circle.  Therefore 

|  b(o)dT(o)  -  f  ii(o)  -r~  ds  where  b££,  n)  is  harmonic 
Y  Y  n 

and  t(£  +  in)  is  an  analytic  function,  is  identically  zero  for  any 
closed  loop  y  which  can  be  shrunk  to  a  single  point  without  crossing 
any  pole  of  b(£*  n)  and  (£  +  in)  (Proof  at  the  end  of  this  section). 

The  unit  circle  can  be  shrunk  to  the  origin  without  passing  through  any 
pole  of  b((,  n)  or  t((  +in)  and  is  one  such  closed  path.  Therefore 


log  (Wfe1^)  -  Wc 
i2ir 


Wc  ' 


we  have 


fl  .  .  i2irt,  .  ,  i2irt.  .  tl  ,  .  i2irt.  3b  . 

/  b(e  )d  log (C  -  •  )  •  1/  log((  -  •  )  -r-  ds. 

t-0  t-0  ** 

Zn  a  similar  manner,  we  have  ' 

1 1  i2irt..  .  i2wt.  fl  .  .  i2irt.  i2irt. 

J  e  db(e  )  -  -  J  b(e  )d(e  ) 

t-0  t-0 

.  -  J1  ,12rt|£d. 

Jt.o  8" 

Consequently  if  the  periodic  function  a(t)  has  its  real 
part  of  bounded  variation  and  its  imaginary  part  continuous  formula 

(10)  becomes 

*<C)  -  /lei2irtd  o(t)  -  Z(C)  -  /1ei2lTtd [Re  (o (t) )  +  2ir|£  t] 

0  0 

-  exp  J^logtC  -  el2irt) d [Re (a (t ) )  +  (1  +  2ir||)  t]  (11) 

where  b(c,  n)  is  a  harmonic  function  defined  on  the  disc  |c|  3  1  and 
equal  to  Im  (a(c))  for  |c|  -  1. 

Putting  ^  -  1/  c  »  W(iJ>)  -  W  -  z(l/*)  -  zo  we  have  a  formula 

o 

corresponding  to  (10)  for  the  internal  problem. 

.  log(W0|0  -  w  )  «  -  J1  log(i  -  ii27rt)d[Y(t)  +  t]  '  (12) 

0  * 

where  y(t)  is  a  complex  valued  periodic  function  of  t  with  bounded 
variation. 

Notes  to  Section  2 

Let  an  integral  Z  be  defined  by 

Z  -  f  b(o)dT(c)  -  [b(o)T(o)]  -  f  T(c)db(o) 

y  Y  Y 

where  y  is  a  closed  path  containing  no  pole  of  either  b(o)  or  t(o),  b(o) 
being  a  real-valued  harmonic  function  of  (£,  n)  ■  (Re(o),  Zm(a))  and 
t(o)  being  a  (complex  valued)  analytic  function  of  a  inside  the  curve  y. 
The  integrals  are  defined  when  one  of  the  functions  is  continuous 
and  the  other  is  of  bounded  variation  on  y. 


We  have 


3t  (g,n)  3t. (g,n) 

I  -  i  bU,  n)  (  -%r-  ♦  i  %.■  ■■)  as. 

Y 

By  the  analyticity  of  t (a)  “  Tr  +  we  can  write 
3t  (5,n)  3T.U*n) 

l  -  if  b u,n>  (-Ejj —  ♦  - 5J— i ,d*  • 


The  real  and  imaginary  part  of  the  above  integral  is  of  the 

,3  C 

form  B — ds  where  B  and  C  are  harmonic  and  n  is  the  normal  for  the 
■lY3n 

curve  y.  Noting  that  Green's  identity  is 


/  (BVC  -  CVB).n<Ss  ■  /  (BVaC  -  CV*B)da  “  0, 

3A  ~  A 

where  A  is  the  part  of  the  plane  <g ,  n)  enclosed  by  the  curve  y,  we  have 

$  b(g,r,)dT(g  +  in)  -  if  t(g  +  in)|^  ds  »  0  (13) 

Y  Y 

for  any  closed  path  y  not  containing  any  pole  of  b(£,  n)  or  x(o),  and  8 
is  the  arc  length  of  the  curve  y . 

As  an  example  we  choose  b(£,  n)  such  that  it  vanishes  on  the 
unit  circle.  The  familiar  real  valued  Green's  function  becomes. 


b(o)  «  log 


1  -  1/a 
1  -  a 


+  log 


o  -  a 
o  -  1/a 


Its  poles  are  at  a  (|a|  <  1}  and  1/a  which  is  outside  the  unit  circle. 
Choosing  the  curve  y  encircling  the  point  a,  then  along  a  radial  line, 
the  unit  circle  and  back  to  the  small  circle  around  a  we  have 


f  b(o)dT(o) 
Y. 


f  b(o)dt (o)  ■  i^ 


T(C)  f  dS 


i$  t(o)  ds. 

u 


The  second  integral  is  zero  as  b(o)  is  identically  zero 
on  u).  The  first  integral  vanishes  as  the  radius  y  reduces  to  zero 


as  t(o)  is  analytic  at  the  point  a.  The  left  hand  side  is  thus  zero. 
The  right  hand  side  gives 


i2n  (a)  ■  i^  t  (o) 
u 


1  ,  c  -  a. 

cos  {arg  - ) 


o-a 


1  .  o 

cos  (arg 


<o-  3) 
a 


1 

a) 


ds 


*  f  t(o) 

(i) 


(o-a)  (o-  i^) 
a 


I)  — 

_  o 
a 


which  reduces  to  a  standard  Poisson  formula  upon  simplification 
(see  [13]  and  its  notes) . 


3.  DETERMINATION  OF  THE  FUNCTION  q(t) 

Let  us  consider  initially  the  problem  of  determining  a(t) 
of  equation  (10)  for  a  given  function  z(C)  of  the  form 

n  r  n 

z(c)  -  z_  +  n  (c  -  0  rk,  E  F  -  1  ,  (14) 

k=l  K  k«l  k 

where  rk's  may  be  complex  and  each  0k  is  inside -the  unit  circle  (|0k|<l) 
unless  its  corresponding  Tk  has  a  positive  integer  value  (1,  2,  3,...). 

In  this  way,  the  function  z(£)  is  defined  unambiguously  outside  the  unit 
circle  |s|  =  1.  Its  derivative  is 


dz 

dC 


2n-l 

n  u-0k 

K«1 


\  2n-l 

.  S  \  «  o, 

k-1 


(15) 


with  all  the  0k  assumed  to  be  contained  in  the  unit  disk  unless  its 
corresponding  £k  is  equal  to  zero.  The  (n-1)  values  of  0k  for  k  =  n  +  1, 
2n  -  1,  are  the  roots  of  a  polynominal  resulting  from  the  differentiation 
of  z(?).  When  two  conditions  of  (14)  and  (15)  are  satisfied  we  have  all 
our  requirements  (1)  satisfied  with  z^  standing  for 

fl  i2irt  . 

z  =  I  e  da (t) .  (16 

c  0 

Now  let  the  function  a(t)  be  the  value  of  a  complex  function 
n(C)  on  the  unit  circle  c  *  el2lTt.  The  problem  is  to  determine  a 
function  n(S)  such  that  equation  (10)  is  satisfied  on  the  unit 
circle  of  the  t  plane. 


jf 


s  /  /  ssy, 


Since  the  contour  integral 


2n-l  . 

j  (  n  (o  -  0 .  )  'nt  “  V  do 

ui  k*l 


has  the  purely  imaginary  value  of  2ni  we  can  define  a  single-vali 
analytic  function  n(C) 

1  f  ,c  2a’1  (Ak  -  V  1 

n(?)  “  2^r  W  1  11  (0  "  0k}  K  )  do  -  log  c  V  , 

o=l  k=l 


where  the  integral  path  is  outside  the  unit  disc.  On  the  unit  ci 
we  have 


1  ,ei2l,t  2n_1  (Ak  "  rk ] 

</  ,  <  n,  (°-v 

o=l  k=l 


Substitution  of  this  function  a(t)  into  equation  (10) 


noting  that 


(\-V  d 


n  <0  -  0  )  K  k  =  §-  (log  n  (o  -  0  )  k) , 


does  give  equation  (14).  Also  substitution  of  this  function  o(' 
into  equation  (8)  gives  equation  (15).  Therefore  the  functio 
a(t)  is  the  required  function. 

As  an  application  of  the  above  calculations  we  can  wr 
down  the  expression  for  (?  -  0)  where  ( 0 |  <  1  as 


C  -  0  *  exp  I/1  log  (?  -  ei2lTt)  yy  d  flog  (el2irt  -  0))j. 


The  familiar  transformation  from  a  circle  to  an  ellif 
or  a  Joukowski  aerofoil  therefore  can  be  written  as 


z (?)  -  =  U  -  m>  (;  -  n)  (?  -  ^y1)"1,  |m| ,  |n|  S  1 


Ul  .  ,  i2nt,  1  Jr.  (el2,It-m 

!°g  (?  -  e  )  —  d [log  - ryy 


If  we  employ  the  method  of  deriving  equation  (11)  then  the 
quantities  inside  the  square  brackets  in  the  last  two  integrals  can  be 
replaced  by  two  purely  real  ones.  However  it  is  simpler  in  practice 
to  use  finite  factorization  rather  than  the  formulae  just  derived. 

The  next  case  we  consider  is  very  practical.  It  is  the 
problem  of  determining  the  transformation  (10)  when  the  shape  of 
a  non-self-intersecting,  continuous,  closed  curve  C  in  the  z-plane  is 
given  in  parametric  form 


w(u)eC,  w(0)  =  w(l) ,  OSuSl 


(19) 


(An  example  is  the  ellipse  w  =  cos(2irua)  +  2isin(2irua)  given  in  the 
z-plane) . 

The  problem  is  to  determine  a  function  z(?)  such  that 


zU)  =  C  +  E  an? 

n=0 


and 

w(u)  =  kz(£)  for  C  =  e'*'2lTt, 

where  k  is  a  positive  scaling  constant  and  t  is  a  monotonic 
continuous  function  of  u  with 


(20a) 


(20b) 


t(0)  =  0,  t(l)  =  1.  (20c) 

For  the  moment  we  assume  that  such  a  function  t(u)  and 
the  function  z(^)  has  been  determined. 

Take  a  point  zc  interior  to  the  image  of  the  unit  circle  | c |  =  1 

and  examine  the  function  (z(?)  -  z  )/?.  This  function  is  single-valued 

c 

and  analytic  for  all  values  of  5  outside  the  unit  circle.  The  image 
of  any  continuous  curve  external  to  |?|  =  1  in  the  C  plane  can  be 
continuously  shrunk  to  a  single  point  +1  in  the  (z(C)  -  zc>/C  plane.  Therefore 
a  single-valued  function  log  ((z(£)  -  zc>/?)  can  be  defined  for  all  points 
|c|21  and  has  a  Laurent  expansion 

2,0  -  '=,  "  ;-n 

log  (  }  )  *  E  ®n 

^  n=0 


with  b  equal  to  zero  due  to  the  condition  (20a) . 


By  taking  a  contour  integral  around  the  unit  circle  |t|  >1 

we  have 

J1  n  (log(zU)  -  z  )-i2irt)e*2lnntdt{u)»0  for  m  *  0,1,2, 3..  (21 

u*o  c 

Especially,  for  m  ■  0 

f\  (log  (w  (u)  -  w  )-i2irt)dt  (u)  -  log  k. 

UBU  c 

Thus 

/J  log|w(u)  -  wc|  dt(u)  *  log  k  (22 

for  all  points  wc  interior  to  the  curve  C. 

The  function  t(u)  determined  by  this  equation  (21)  is  the 
required  function  if  it  is  monotonic.  It  is  not  difficult  to  see  that 
dt/du  is  the  charge  distribution  (with  respect  to  u)  on  the  boundary 
C  of  a  two  dimensional  conductor  in  an  external  potential  problem. 

Once  the  function  t(u)  is  solved  from  (22),  the  required 

function  a(t)  «  log  (z(el2lrt)-  z  )  is  determined.  There  is  only 

c 

one  function  t(u)  satisfying  equations  (20).  Its  uniqueness  is 
proved  by  showing  that  any  function  z(£)  satisfying  all  conditions 
( 1 )  and 

|ztt)|  -  1  for  |c|  =  1 

must  be  identically  equal  to  £.  The  proof  uses  Laurent  series 
expansion  and  the  Maximum  Modulus  theorem. 
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THE  RESTRICTION  POSED  BY  EQUATION  (7) 

The  restriction  (7)  has  to  be  observed  for  the  proposed 
transformation  (10)  to  work.  The  physical  meaning  of  this  restriction 
is  now  examined.  Equation  (11)  has  shown  that  for  any  continuous 
a(t)  we  can  find  another  real  function  [Re(a(t))  +  (1  +  2n9b/3n)t] 
of  equivalent  effect  therefore  we  shall  content  ourselves  with  the 
case  of  only  real  a(t)  in  this  section. 

Consider  the  real  part  of  the  last  integral  on  the  right 
hand  side  of  equation  (10) .  It  is  the  potential  generated  by  a 
charge  distribution  proportional  to  da/dt  on  the  periphery  of  the  unit 
circle.  The  restriction  (7)  then  requires  that  its  corresponding 
electric  field  should  not  vanish  anywhere  outside  the  unit  circle. 

This  vanishing  can  precipitate  the  formation  of  kinked  potential 
lines  (such  as  those  apple  shaped  curves  associated  with  log  z  * 

(1  +  k)  log  c  -  k  log  (c  -  a))  which  leads  to  a  finite  area  in  the 
z  plane  being  the  image  of  a  line  segment  in  the  C  plane.  In  other 
words  the  area  of  uniform  potential  in  the  z  plane  is  more  than  the 
area  interior  to  the  curve  C#  hence  there  may  exist  some  area  of 
uniform  potential  adjacent  to  some  two-dimensional  conductors  which 
is  not  totally  enclosed  by  the  conductors.  There  are  analogous 
results  for  three  dimensional  problems .  One  advantageous  use  of  such 
results  is  with  the  operation  of  Van  der  Graaff  static  generators. 

There  is  more  than  one  way  to  satisfy  the  condition  (7) 

The  easiest  way  to  satisfy  this  condition  is  to  have  (a(t)  +  t] 
monotonically  increasing. 

The  existence  of  those  "regions  of  degeneration”  has 
diverse  physical  meanings  besides  that  already  cited:  In  problems 
of  potential  flows  we  have  regions  of  separation,  in  wave  scattering 
problems  we  require  points  of  complex  mass  attached  to  the  boundary 
of  the  basic  circular  hole  to  generate  all  patterns  of  far-field 
scattering. 

5.  RELATIONSHIP  TO  THE  SCHWARZ -CHRISTOFFEL  TRANSFORMATION 

Remember  that  a  Schwarz-Christoffel  transformation  is  of  the 
form  given  by  equation  (4).  This  makes  the  point  |  £  |  =  <*  a  singular 
point.  We  now  consider  a  modified  formula 


4-  1  ®  ^ 

aT  “  Ta  n  (c  -  0k)  ,  Yk'»  «U  |0kl  ■  1,  2  Yk  “  1  (23) 

^  k»l  k*l 


Zt  is  obvious  that  s(()  describes  a  polygonal  curve  when  c  describes 
the  unit  circle  just  as  with  the  Schwarz-Christoffel  transformation. 
However  the  transformation  (23)  can  satisfy  condition  (lc)  whereas 
the  other  can  not. 

By  writing 


„  u  n  i2sk/m.m 

(  ■  lim  n  (c  -  e  ) 

m-H»  k»l 


we  can  readily  put  equation  (23)  into  the  form  (15)  which  has  been 
proved  to  be  a  particular  result  of  the  formula  (10) . 

Using  a  to  denote  1/t  we  can  mite  (23)  as  dz/dc  ■ 


jJIj  (l/o  -  0fc)  ,  which  gives 


|  II  t  (^  -  0R)  (o  -  0k)l  *  }  ■  n  (o  *  &k>  k  on  the  circle  |o|  »  1 

l  k«l  k»l 


This  last  expression  is  very  similar  to  the  Schwarz-Christoffel 
formula  (4) .  we  can  therefore  consider  the  latter  transformation  as 
a  composite  one  from  the  o  ■  1/c  plane  to  the  final  S  plane  with  the 
ratio  dS/dz  real  on  the  circle  |o|  *  1.  This  consideration  explains 
its  peculiar  mapping  from  the  interior  to  the  exterior  of  the  unit  circle 
in  the  a-plane  as  n  tends  to  infinity. 

Just  as  with  the  Schwarz-Christoffel  transformation  there 
is  a  linear  form  for  formula  (10) .  This  is  obtained  by  applying  the 
inverse  Mobius  transformation  t  *  i (1  +  C)/(l  -  O  to  the  (  plane. 

The  combined  transformation  z(c(t))  is  then  from  the  half-plane 
Zm(T)  <0  onto  the  exterior  of  a  closed  continuous  curve  C. 


6.  APPLICATION  OF  FORMULA  (6)  TO  THE  GENERATION  OF  AEROFOIL  SECTIONS 

In  generating  aerofoil  sections  it  is  preferable  to  employ 
formula  (6).  The  method  is  to  have  one  zero  of  order  (1  -  b)  for 


[ 


dz/dt  on  the  unit  circle  | ( |  ■  1  while  keeping  all  its  other  singular 
points  inside  the  unit  circle.  The  real  value  b  is  the  internal  angle 
of  the  corner  (trailing  edge)  of  the  aerofoil  section. 

The  Joukowski  sections  are  thus  the  case  with 
*(C)  -  zc  -  ((  -  a)(C  -  n>  tt  -  21S)"1  ,  (26) 


dz 

dc 


m+n  jin-n  m+n 

r  ~  2  *  1  2  )  r  ~  2 


(5 


m+n.  a 
2  ' 


m^n 
2  ) 


(27] 


and  the  complex  constants  m,n  are  chosen  such  that  one  of  the  values 
t  i— is  on  the  unit  circle  and  the  other  is  inside.  The  pole  of 
dz/dc  is  the  center  of  these  two  points  and  is  therefore  put  inside 
the  unit  circle. 


The  Karman-Trefftz  sections  are  given  by 


q(o) 


( 


P  +  1 

o-l 


)  +  1 


(° +  V  - 

o  -  r 


1  <  V  S  2 


(28: 


and 


da 

do 


_ 4v* 

(o  +  1)  (o  -  1) 


( 


0  +  l.v 

o  -  l' 


0  +  l.v 
0  -  lJ 


a 


(29: 


which  require  cuts  between  q  ■  -v  ,  q  »  +v  and  o  ■  -1,  o  ■  +1  in  the 
q  and  o  planes  respectively  (see  [14] ,  page  64) . 


have 


Employing  the  intermediate  transformations  w(q)  and  u(o)  we 


( 


o  +  l.v 
o  -  l' 


This  shows  clearly  that  there  is  no  singularity  of  dq/do  outside  the 
circle  |o|  61  save  for  a  possible  singularity  at  |o|  *»  .  This 
latter  point  is  shown  by  to  be  not  a  singular  one.  Two  new 

variables  z(q)  and  ((o)  are  then  introduced  as  the  same  linear  function 
of  q  and  o  respectively.  The  composite  transformation  z(£)  then  gives 
a  Karman-Trefftz  section  as  the  image  of  a  unit  circle  |c|  -  1. 


Note  that  the  image  a(c)  of  ](|  *  1  in  the  o  plane  must  not  cross  the 
cut  between  a  *  -1  and  a  ■  +1.  By  examining  the  intermediate  variables 
u(a)  and  w(q)  we  see  that  the  last  condition  automatically  ensures  that 
no  singularity  for  z(0  can  exist  outside  the  disc  |c|  S  1.  • 

Following  the  previous  two  methods  of  generating  aerofoil 
sections  we  can  have  quite  a  range  of  transformations  in  the  form 

al  °2  “k  k 

2(0  -  *c  »  U  +  fflj)  (c  +  m2)  ...(c  +  s^)  ,  E  ■  1,  (3 

where  cr's  are  all  real.  For  »(0  to  be  defined  uniquely  every  m^  must 
be  contained  inside  the  unit  disc  | ( |  21  unless  its  corresponding  power 
index  is  a  positive  integer  {1 ,  2 ,  3  . . . ) . 

The  derivative  dz/d(  is  given  by 


dz  1 

a7  *  <*<C)  -  *  >  [— — —  + 


+  . . .  + 


Ar  l.  ,  _  ’  ,  .  _  T  •••  '  .  ,  _  I 

“C  c  C  +  C  ▼  ®2  C  + 

and  all  singularities  (zeros  and  poles)  of  dz/dc  are  contained  within 
the  disc  |c|  2  1.  This  last  condition  in  turn  requires  that  |m^|  is 
inside  the  unit  disc  |(|  21  unless  its  corresponding  power  index 
is  equal  to  one.  (The  Joukowski  aerofoil  given  by  (6.1)  has  one 
zero  of  order  one  of  z(0  outside  the  unit  circle). 

For  examples,  a  few  aerofoil  shapes  have  been  produced 
using  the  following  formulae 


z(c)  -  se  -  (C  +  D  U  -  0.1  -  0.2i)  U  +  0.1  -  0.2i) 


2(0  -  ze  -  (t  ♦  1)  U  -  0.1  -  0.21)  U  +  O.li) 


1.9  -0*45  -0*45 

2(0  -  -  U  +  1)  U)  U  ♦  0.1)  U  , 


or  we  can  even  use  complex  powers  (in  this  case  we  require  that 

|m. I  21  if  (a.  -  1)  is  real  and  non-zero, |m, I  <  1  if  a.  is 

i  i  •  •  ‘  i '  i 

complex)  to  have 


*U>  -  *  «  U  ♦  l)1,9(c)0,1"0,3i(c  +  0.1)‘1+0*3i. 


Not  all  choices  of  m^  and  give  acceptable  aerofoil  shapes. 
For  the  above  four  functions  z(()  a  computer  program  has  been  used  to 
ensure  that  all  singularities  of  dz/dc  are  contained  in  the  unit  disc 
|c|  S  1  except  for  one  at  (  ■  1  to  generate  a  finite  acute  angle>  at 
the  trailing  edge.  The  resulting  aerofoil  sections  have  been  plotted 
by  a  computer  and  the  effect  of  m^'s  and  a^’s  on  thickness,  camber 
etc.,  was  examined  (The  computer  programs  are  written  by  Mr.  C.  A. 
Martin,  the  author  acknowledges  his  help  to  substantiate  the  proposed 
method) . 

7.  CONCLUSIONS 

The  transformation  (10)  subjected  to  condition  (7)  is 
a  continuous  form  of  that  given  by  equation  (23) .  Our  representation 
for  x(C)  in  equation  (6)  or  (10)  has  some  similarities  to  Wierstrass' 
theorem  on  entire  functions  (see  [10] ,  page  195) . 

The  shapes  not  covered  by  the  transformations  (6)  or 
(10)  due  to  the  restriction  (7)  may  have  regions  of  finite  area 
which  are  not  in  the  range  of  the  conformal  mapping.  For  these  cases 
it  is  mathematically  satisfactory  to  modify  the  boundary  to  exclude 
those  regions. 

It  is  noted  that  the  conformal  mapping  here  has  some 
relationship  to  the  more  general  problem  of  the  existence  of  a  conformal 
mapping  from  the  interior  of  a  unit  circle  to  an  arbitrary  simply 
connected  region  (see  [10] ,  page  230) .  The  latter  problem  is  associated 
with  such  famous  mathematicians  as  Riemann  and  Wierstrass  (see  [15] , 
page  186) .  The  difference  between  the  mapping  considered  here  and 
that  of  the  general  problem  is  our  requirement  (lc) .  It  ought  to  be 
mentioned  here  that  the  numerical  method  which  combines  numerical 
solution  to  the  external  potential  field  for  a  given  conductor  shape 
with  a  Fourier  series  analysis  to  determine  the  series  expansion  (5) 
for  the  corresponding  conformal  mapping  can  be  considered  a  logical 
outflow  of  Riemann 's  idea  for  the  existential  problem. 
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